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Abstract 



Extremal dynamics represents a path to self-organized criti- 
cality in which the order parameter is tuned to a value of zero. 
The order parameter is associated with a phase transition to an 
absorbing state. Given a process that exhibits a phase transition 
to an absorbing state, we define an "extremal absorbing" process, 
providing the link to the associated extremal (nonabsorbing) pro- 
cess. Stationary properties of the latter correspond to those at 
the absorbing-state phase transition in the former. Studying the 
absorbing version of an extremal dynamics model allows to deter- 
mine certain critical exponents that are not otherwise accessible. 
In the case of the Bak-Sneppen (BS) model, the absorbing version 
is closely related to the "/-avalanche" introduced by Paczuski, 
Maslov and Bak [Phys. Rev. E 53, 414 (1996)], or, in spread- 
ing simulations to the "BS branching process" also studied by 
these authors. The corresponding nonextremal process belongs to 
the directed percolation universality class. We revisit the absorb- 
ing BS model, obtaining refined estimates for the threshold and 
critical exponents in one dimension. We also study an extremal 
version of the usual contact process, using mean- field theory and 
simulation. The extremal condition slows the spread of activity 
and modifies the critical behavior radically, defining an "extremal 
directed percolation" universality class of absorbing-state phase 
transitions. Asymmetric updating is a relevant perturbation for 
this class, even though it is irrelevant for the corresponding nonex- 
tremal class. 
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I. INTRODUCTION 



Extremal dynamics has been employed extensively in modelling far from 
equilibrium systems such as biological evolution [1] and driven interfaces [2,3]. 
Although processes with extremal dynamics do not have a phase transition 
(there is no control parameter) they exhibit scaling properties reminiscent of 
those observed at continuous phase transitions [4,5]. Indeed, it was suggested 
some time ago that the appearance of "self-organized" scaling properties in 
extremal dynamics and in sandpile models corresponds to forcing the order 
parameter (associated with an underlying phase transition) to zero from above 
[6]. The purpose of this work is to explore the connection between these scal- 
ing properties and those observed at a phase transition to an absorbing state. 
The connection between extremal dynamics and directed percolation (DP), 
the prime example of an absorbing-state phase transition, was first suggested 
by Paczuski, Maslov and Bak [4] and investigated in detail by these authors 
in the context of the Bak-Sneppen (BS) model and related processes [5]. (The 
latter work, as well as Ref . [7] , clearly demonstrated that the critical exponents 
of the BS model are not those of DP.) Sornette and Dornic [8] and Grassberger 
and Zhang [9] have shown how a variant of directed percolation may be trans- 
formed via extremal dynamics to display SOC. These studies indicate that 
self-organized criticality (SOC) [10] under extremal dynamics arises because 
the system is driven to a critical point associated with a phase transition to 
an absorbing state [8], as is also the case for sandpiles [11]. 

In the present work we are particularly interested in the modifications 
needed to transform a (non-extremal, non-SOC) model having an absorbing 
state to one exhibiting SOC under extremal dynamics. We develop a gen- 
eral scheme relating the two classes of models via an intermediate, "extremal- 
absorbing" process whose absorbing-state critical point corresponds exactly to 
the critical behavior observed in the corresponding SOC model. Two examples 
(the BS model and an extremal contact process) are studied in detail, yield- 
ing refined estimates for critical properties, and evidence of a new universality 
class associated with absorbing phase transitions under extremal dynamics. 

The prime example of extremal dynamics is the BS model [1,12], proposed 
to explain mass extinctions observed in the fossil record. While its application 
in the evolutionary context is debated [13], it remains an intriguing and in- 
completely understood example of scaling behavior far from equilibrium. The 
contact process (CP) [14] is the most familiar example of a Markov process 
exhibiting a phase transition to an absorbing state. We focus on the absorbing 
version of the BS model, and the extremal version of the CP, to illustrate the 
relations between extremal dynamics and absorbing phase transitions. Our 
analysis of the spread of activity leads to new or refined values for the expo- 
nents 5, i], j3' and z sp , and for the critical threshold of the BS model. (For 
extremal dynamics, the avalanche exponent is r = 1 + 5.) Finite-size scaling 
analysis of stationary properties at the critical point yields estimates of the 
exponent ratios P/v± and v\\/v±- 
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Studies of modified BS models have shown that scaling properties are insen- 
sitive to changes that preserve its basic symmetries (that is, invariance under 
translation and reflection) [15-17], pointing to the existence of a BS universal- 
ity class. Nonextremal models that exhibit a phase transition to an absorbing 
state, and that possess these same symmetries, and no additional ones, belong 
generically to the directed percolation (DP) universality class [18,19]. Here we 
show that such models fall in a new "extremal-DP" universality class when 
modified to follow extremal dynamics. The absorbing phase transition corre- 
sponding to the scaling behavior of the BS and other extremal models belongs 
to the extremal-DP class, not that of ordinary directed percolation. In other 
words, extremal dynamics is a relevant perturbation for absorbing-state tran- 
sition, just as was shown by Sneppen in the context of interface depinning 
[2]- 

The balance of this paper is organized as follows. Section II presents a 
general scheme linking ordinary absorbing-state phase transitions and extremal 
dynamics via an intermediary "extremal-absorbing" model. In Sec. Ill we 
describe our simulation method and results. Our findings regarding scaling 
and universality are discussed in Sec. IV. Mean-field analyses are presented in 
the Appendix. 

II. ABSORBING STATE MODELS AND EXTREMAL DYNAMICS 

In this section we examine how a stochastic model with an absorbing- 
state phase transition may be transformed to exhibit scale invariance under 
extremal dynamics. We begin, for generality, by defining a rather abstract 
scheme, and then discuss specific examples. A large class of models exhibit- 
ing an absorbing-state phase transition may be formulated as follows [20-23] . 
Consider a stochastic process S defined on a connected graph Q of N sites. (Q 
consists of a set of sites with links between certain pairs of sites. Typical ex- 
amples are a ring of N sites, and the d- dimensional hypercubic lattice Z d , with 
links between nearest neighbors.) The state a{i) of site i is or 1, the latter 
value denoting an active site, the former an inactive one. For each site % in Q 
we define a neighborhood v(i) C Q, or, more generally, a set of neighborhoods 
vi(i), v 2 {i),...,v n (i). 

The dynamics of S proceeds in steps. Each step involves choosing an 
active site i (the central site for this step), at random, and changing the states 
of the sites in v(i) according to a certain rule /. In case there are two or more 
neighborhoods, one of them, v r {i) say, is chosen at random from the collection, 
with probability p r , and a rule f r is applied to the site or sites in v r . In general 
/ (or f r ) is a probabilistic rule. At each step the number of active sites may 
change, and if at any moment there are no active sites (cr(i) — 0, Mi e Q), the 
process has fallen into an absorbing state and there is no further evolution. 
Otherwise the dynamics proceeds to the next step. 

Using a n to denote the entire set of activity variables cr(l), cr(2), cr(N) 
at step n, the dynamics generates a sequence <7i, <r 2 , ... starting from the initial 
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configuration a . It is frequently of interest to associate a continuous time vari- 
able t with the process. This is usually done by associating a time increment 
At = 1/N a with each step, where N a is the number of active sites just before 
the step is realized. We define the order parameter as pit) = Prob[a t (i) = 1], 
i.e., the fraction of active sites at time t. (The event space here is the set of all 
realizations of the process up to time t, starting from a given initial probability 
distribution on configuration space.) If the stationary order parameter, defined 
by linn^oo hin^^oo p, vanishes, the process is said to be in the absorbing phase; 
otherwise it is in the active phase. 

A simple model exhibiting a phase transition to an absorbing state is the 
contact process (CP) [14]. Here we consider the one-dimensional version. There 
are three sets v r (i), conveniently denoted as vq(i) = i and v±(i) — % ± 1. The 
associated probabilitites are p and p± = (1 —p )/2. In terms of the usual 
parametrisation [14,23], p = 1/(1 + A), where A > represents the rate of 
spread of activity. (In the "epidemic" interpretation of the CP, active sites 
represent infected organisms, inactive sites susceptibles, and A is the infection 
rate.) The updating rules are: fo = 0, f± = 1. In other words, an active 
site has a probability per unit time of 1/(1 + A) to become inactive, while an 
inactive site j becomes active at rate \n a (j)/ [2(1 + A)], where n a (j) is the 
number of active neighbors of site j. As is well known, the one-dimensional 
CP exhibits a continuous phase transition between an absorbing phase and an 
active one at A c ~ 3.29785 [20,22,23]. 

It is convenient to associate the control parameter with the updating rule 
/ rather than with the probabilities p r . We therefore reformulate the CP as 
follows. With the sets v (i) and v±(i) defined above, we take p = 1/2 and 
p± = 1/4, and define q — X/(l + A). The updating functions are: 



It is easy to verify that the transition rates satisfy w(0 — > l)/w(l — > 0) = 
n a \/2, just as in the original formulation. The critical value q c ~ 0.76733. 

The following three-site contact process (CP3) will play an important role 
in our analysis [24]. For each site we define the set v(i) = {i — l,i,i + l} (the 
central site and its nearest neighbors). The updating function / takes values of 
1 and with probabilities q and 1— q respectively, independently at each of the 
three sites in v(i). (In Ref. [24] this is called 'model 3'.) Simulations of the CP3 
show that it exhibits a continuous phase transition at q = q c ~ 0.63523(3). 
Further interacting particle systems, such as the pair contact process [25] and 
the diffusive CP [26], can be accomodated within the scheme set out above. 




(1) 



('w.p.' denotes "with probability"), and 
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We shall assume that the process S is defined so that the control parameters 
(such as q) are associated with the updating rule /. Each time a site is updated, 
the value of / may be determined by comparing a random number x with the 
parameter in question. (This is of course the usual procedure in simulations.) 
In the CP3, for example, we take / = 1 if x < q, and zero otherwise, where x is 
uniformly distributed on the interval [0,1]. Call the random number associated 
with the most recent updating of site i, Xi, so that a(i) = Q(q — Xi) with the 
unit step function. (The initial values of the Xi are assigned according to the 
state variables cr(i). For example, if all sites are initially active, we draw the 
initial X; L from the distribution uniform on [0, q\.) For the CP, Eq. (1) requires 
that we update the central site % with a number chosen uniformly from [0,1]. 
According to Eq. (2), the same applies when updating an inactive neighbor 
(i ± 1), but when updating an active neighbor, the random number x is drawn 
from the interval [0, q], since an active neighbor remains active. 

Summarizing, we have shown how a particle system S may be formulated 
using a set of random variables Xi, such that site % is active if Xi is smaller than 
a certain parameter q. S suffers an absorbing-state phase transition at q — q c . 
We now define two related processes, Sea and Se- The former, extremal 
absorbing process, is obtained by modifying how the central site is selected. 
Instead of choosing it at random from among the currently active sites, it is 
taken to be the active site having the smallest Xi. As in the original process 
S, if there are no active sites (i.e., Xj > q, Vj G G), the process has reached 
an absorbing configuration and the evolution ceases. Thus Sea possesses an 
absorbing state, and since the original process S exhibits a phase transition 
between an active and an absorbing phase as the control parameter q is varied, 
we expect Sea to as well, at some value q C} E- (The reason is that the relative 
likelihood of generating and destroying active sites varies with q, just as in 
S.) Mean-field theory (see Appendix) yields q C) E = q c - F)ue to the different 
correlations generated under extremal dynamics, however, the critical value 
q c> E of Sea is in general different from q c . 

We define the extremal process Se by relaxing the condition in Sea, that the 
extremal site % must be active (i.e., have Xi < q) for the dynamics to proceed. 
If S is the original contact process, then Se is a process in which either the 
minimal site or one of its nearest neighbors is updated at each step. If S is the 
CP3, S E is the familiar Bak-Sneppen model. Note that S E has no absorbing 
state, hence no phase transition to such a state. Its stationary properties are 
nevertheless intimately connected with the critical-point properties of Sea, as 
we now explain. 

Of particular interest is the stationary probability density p(x) of site vari- 
ables under extremal dynamics. As is well known, p(x) = C6(x-g c ,B)0(l-i) 
in the Bak-Sneppen model, in the infinite-size limit. (C = 1/(1 — q c ,E) is the 
normalization factor.) We expect p(x) to exhibit a step-function singularity 
in any extremal model Se [15]. This feature is in fact already present in the 
original absorbing-state model S at its critical point, because at the critical 
point q = q c , the stationary density of active sites (having Xj < q) tends to 
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zero as the system size N goes to infinity. The distribution on the "allowed" 
region x > q is uniform, since the Xj are drawn from a uniform distribution. 
Thus p(x) jumps from zero to a finite value at x — q c . In the supercritical 
regime (q > q c ), p(x) is equal to a constant p\ for x < q, such that qpi = p 
(the order parameter), and takes a different constant value, P2, on the interval 
[q, 1]. Once again, the stationary density is discontinuous at x — q. 

What holds for S also holds qualitatively for Sea- The critical value q C} E 
may, as noted, differ from q c , but since Sea exhibits an absorbing-state phase 
transition, its stationary distribution p(x) also has a step-function singularity. 
Just at q — q C) E, the order parameter p = 0, but as iV — > oo the survival 
time of the process tends to infinity. This means that the process can survive 
indefinitely, with the choice of the central site restricted to the set having 
x < q c ,E- The presence of active sites in the range g Cj £ < x < q is then 
irrelevant, since a site with x < q Ct E is always available. Thus for q > q c ^E, 
the distribution p(x) exhibits a step-function singularity at q c> E- Extremal 
dynamics effectively "pins" the singularity at q Cj E- The foregoing remarks on 
Sea obtain in the infinite-size limit; for finite N there is a nonzero probability 
(for q < 1), that all sites have x > q so that the system eventually becomes 
trapped in the absorbing state. (The mean lifetime, however, is expected to 
grow exponentially with N, for q > q Cj E-) 

In summary, if S exhibits an absorbing-state phase transition, then Sea 
should as well, although not necessarily at the same value of q. The distribution 
p(x) possesses a step function singularity in both cases. Near the critical 
point (q > q c ) of the original process S we expect p ~ (q — qc) 13 , with (5 
the critical exponent associated with the order parameter. Below the upper 
critical dimension (d c = 4 for DP [18,19]), (5 < 1. In the supercritical regime 
of Sea, on the other hand, p = f^p(x)dx oc (q — q c ,E) since p(x) jumps from 
zero to a finite value at x — q c ,E- Thus the order parameter exponent (3 is 
unity in Sea- This is but one example showing that the extremal condition 
changes the value of a critical exponent. Another example, established some 
time ago by Paczuski and coworkers [5,27], is that the spreading exponent 77 
(defined Sec. III.C) is generically zero for extremal dynamics, whereas rj > 
for directed percolation. Further evidence that extremal dynamics modifies 
critical exponent values will be given below. 

In Se, the central site is always chosen (in the iV — > 00 limit) from the set 
{i : Xi < q c ,E}, just as in Sea at its critical point. Since the order parameter is 
zero in the latter case, we may assert that extremal dynamics effectively tunes 
the order parameter to zero, (p approaches zero from above as N — > 00.) As 
we have seen, the existence of sites with q c E < x < q becomes irrelevant in 
Sea, in the infinite-size limit. In other words, Se is identical to the critical 
process Sea in this limit. (Starting from the same initial configuration, and 
using the same set of random numbers, the same sequence of sites will be 
updated in the two processes.) 

Since the extremal version of the CP3 is the familiar Bak-Sneppen model, 
we shall refer to the CP3ea as the absorbing Bak-Sneppen (ABS) model. Our 
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objective is to characterize the behavior of the ABS model and of the extremal 
and EA versions of the contact process. The ABS model is closely related to 
the f- avalanche process studied in Ref. [5]. An /-avalanche (in the present 
notation, g-avalanche), begins when the minimal site variable x min < q, the 
minimum having been greater than q at the preceding step or steps, and con- 
tinues until the minimum is once again > q. (As q approaches q C) E from below, 
the mean avalanche duration diverges.) The dynamics of the BS model contin- 
ues, regardless of whether a given avalanche has terminated or not. But in the 
ABS model x min > q represents an absorbing state and the dynamics ceases. 
In the BS model, it is common to analyze the properties of g-avalanches in the 
stationary state. It is similarly of interest to study stationary properties of the 
ABS model, attained once the system has relaxed, after an initial transient pe- 
riod. We may also study the mean lifetime of the active state as a function of 
system size. Another approach to studying absorbing-state phase transitions 
consists in following the spread of activity starting from a single active site. 
This spreading phenomenon in the BS model was studied in Refs. [4,5], where 
it is called the BS branching process, and in Ref. [7] under the name of the 
BS(p) model. 

The assertions regarding extremal and extremal-absorbing models are sup- 
ported by mean-field theory (MFT), which is presented in the Appendix. In 
particular, for the ABS model p = §0(x — |), just as in the MFT of the orig- 
inal BS model. In the extremal and the extremal-absorbing contact process, 
the stationary probability densities exhibit (for q > q Cy E = 1/2), two disconti- 
nuities, one at x = k = q 2 /(Sq — 1), the other at x = q > k. [These coalesce 
at q — 1/2. Note that the parameter q continues to influence the form of p in 
the extremal contact process, due to the nature of the updating rule, Eqs. (1) 
and (2).] These predictions are in qualitative agreement with simulation. 

III. SIMULATION RESULTS 

A. Absorbing Bak-Sneppen Model: Simulation Method 

Before discussing our results we insert an observation on simulations of the 
BS model. Since the site with the smallest variable, x m i n , must be identified 
at each step, it becomes important to devise an effective search strategy. An 
efficient general-purpose search algorithm uses a binary tree structure to iden- 
tify x min . One approach [7] utilizes a lattice of 2™ sites. At the first level of 
selection, each site is compared with one of its neighbors and the minimum of 
the pair selected. At the next level the minimum between each neighboring 
pair is selected, and so on, so that at the n-th level the global minimum is 
identified. 

A second binary scheme [28] is formulated as follows. Site is placed at 
the apex of the tree. Site 1 is placed on the level below the apex, to the left 
of if x\ < xo, to the right if x\ > xo- A site i is added to the tree in the 
following way: we go down the tree comparing Xi with the variables x±, Xi-i, 



7 



turning left or right depending on whether Xi is smaller or larger than Xj, until 
we find an empty site. Building the tree in this way, x m i n will occupy the 
leftmost position in the tree. In these schemes, maintaining the tree structure, 
once constructed from the initial set of variables x iy requires a small number 
of operations at each step, and is many times more efficient that a repeated 
global search for the minimim. We find, nonetheless, that a suitably restricted 
search requires less cpu time in the stationary state. 

A special property of the BS model (shared by its absorbing version, and 
by other extremal models), is that the minimal site falls, with a probability 
approaching unity as the system size grows, in the interval [O,^^]. At the 
same time the density of sites with values in this interval approaches zero as 
N — > oo. This suggests maintaining a list of sites having x < q C)E [29]. Then 
the search for x min may be restricted to the list, except for the rare instances 
in which the latter is empty. (For the ABS we must in any case restrict the 
search to sites with x < q.) If the system is large, so that the typical number 
of sites with x < q is not small, it becomes advantageous to introduce a second 
list, of sites having x < q** < q CyE - When this relatively short list is nonempty 
(as is usually the case) the search for x min is restricted to it. In studies of 
the BS model, we obtained the greatest efficiency using q** = 0.54, while 
the criterion for the first list was x < 0.65, that is, slightly below q C:E . (The 
occasional need to perform a global search, in the rare instances when both 
lists are empty, is more than compensated by their reduced sizes when using 
these values.) Compared with the binary tree method, our approach results 
in threefold reduction in CPU time, in the stationary state, for a system of 
1000 sites. (The binary tree approach may prove more efficient for studying 
transients, since initially the lists will not be short, if the Xi values are chosen 
uniformly on [0,1].) 



B. Absorbing Bak-Sneppen Model: Stationary Properties 

Using the simulation method described above, we conducted extensive stud- 
ies of the one-dimensional ABS model. We initialize the system with all sites 
active and allow it to relax until mean properties fluctuate about stationary 
values. The stationary properties are then obtained from temporal averages 
over the set of surviving realizations. Each step corresponds to a time interval 
of At = 1/N a , with N a the number of active sites just before the updating is 
performed. Results for the stationary order parameter (i.e., the density of sites 
with x < q), on lattices confirm that in the supercritical regime (q > q c ,E), the 
order parameter grows linearly with q — q c>E , as anticipated in Sec. II. 

We study the finite-size scaling behavior of the stationary order parameter p 
and of the lifetime r at the critical point, (r is obtained from an exponential fit 
to the survival probability P s (t).) The expected finite-size scaling behaviors at 
the critical point are: p ~ L~P/ U± and r ~ L v \\l v ^. We performed simulations 
at q — 0.66701 and 0.66702, the latter being the preferred literature value for 
the threshold in the one-dimensional BS model, while the former is favored by 
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the results discussed in the following subsection. We studied systems of 1000, 
2000, 4000,. ..,32000, sites in simulations of 2 x 10 7 to 3 x 10 s time steps. 

For L = 4000 - 32000, the results for the order parameter follow a power 
law with p/u± — 0.755(5). The data for smaller system sizes, however, show 
systematic deviations from a pure power law, leading us to seek a correction 
to scaling term; a correction decaying oc L -1 / 2 leads to a good fit. We fit the 
expression 

lnp = -Ai nL __L (3 ) 

to the data for L > 1000, allowing j3/u± and b as adjustable parameters; the 
best-fit values are (3/v±_ = 0.769(7) and b = 3.69(20). The simulation data, 
and the difference from best fit of Eq. (3) are plotted in Fig. 1, showing 
the high quality of fit. The data for the lifetime, using system sizes of 125, 
250, ...8000 yield v\\/v±_ = 2.12(1), with no obvious correction term (see Fig. 
2). We also determined the stationary moment ratio m = p 2 /~p 2 at the critical 
point. The estimates for m decrease slowly with L, and appear, when plotted 
versus L~ a25 , to approach a limiting (L — > oo) value of 1.030(5) (see Fig. 3). 
(For the one-dimensional CP, m = 1.1737 at the critical point.) Essentially 
the same results are obtained regardless of whether we use q CjE = 0.66701 or 
0.66702 in these simulations. 



C. Absorbing Bak-Sneppen Model: Spread of activity 

Scaling properties at an absorbing state phase transition are also reflected 
in the spread of activity from an initially localized region [30]. In spreading 
simulations of the ABS model we start the system with a single site x < q and 
all others above this value. (This is completely equivalent to the BS branching 
process studied in [4,5,7].) At q — q c ,E, the process generates a scale-invariant 
pattern of activity that may be characterized by power-laws for the survival 
probability P s (t), the mean number of active sites n(t) and the mean-square 
distance R 2 (t) = [™(£)] _1 (X]j r 2 (£)). (rj(t) denotes the position of the j — th 
active site at time t. Note that n(t) is taken over all realizations, including 
those that have become trapped in the absorbing state at or before time t.) 
The scaling laws are typically written in the form 

P s ~r s , n~t", i? 2 ~t Zsp , (4) 

relations that have been verified to high precision for various examples [21,23]. 
(We use z sp to denote the spreading exponent, to avoid confusion with the 
dynamical exponent z.) The appearance of power laws is commonly used to 
locate the critical point [23]. 

The spreading exponent 5 is related to the avalanche size exponent r, de- 
fined (in the BS model) via Pd(s) ~ s~ T , where P D (s)ds is the probability 
of an avalanche having a duration between s and s + ds. Thus the survival 
probability P s (t) = J* t °° Po(s)ds, implying r = 1 + 5. 
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We performed spreading simulations of the ABS model at q — 0.66699, 
0.66700, 0.66701, 0.66702 and 0.66703. Each realization was followed up to 
a maximum time of about 2.7 x 10 5 ; the total number of realizations ranged 
from 4 x 10 5 to 1.6 x 10 6 , depending on the value of q. To locate the critical 
point we study the local slope 5(t) = dlnP/dlnt, plotted versus For 
q < q c the local slope is expected to veer downward at large times, and vice- 
versa. (Numerically, S(t) is given by the slope of a least-squares linear fit to 
the data in an interval [t ,20t ], with geometric mean t.) On the basis of 
the local slope data (see Fig. 4) we conclude that q CjE = 0.66701(1). This is 
consistent with previous estimates, which place the threshold at 0.66702(8) [7] 
and 0.66702(3) [5]. (We did not find analyses of the local slopes r)(t) or z sp (t), 
defined analogously to S(t), useful in locating the critical point.) 

Analyzing the data at the critical point, we are unable to obtain good fits 
to P s , n and R 2 using simple power-law expressions. Including a subdominant 
power-law correction in the relations of Eq. (4) greatly improves the quality 
of fit. In particular, the survival probability can be fit quite accurately using 

lnP s ~ -51nt + <f)pt- 1/4 + C (5) 

where C is a constant and the best-fit values are 8 = 0.084(1) and <pp = 0.115. 
The same value for 5 is found using the data for q = 0.66702. (The choice 
of a correction term decaying as t -1 / 4 is motivated by the fact that the local 
slopes 5(t) and z sp (t) are essentially linear when plotted versus t -1 / 4 , as seen 
in the inset of Fig. 4.) In Fig. 5 we plot P s and the ratio of P s to the fitting 
function, Eq. (5); the ratio is seen to be essentially constant for t > 50. The 
mean-square displacement may also be fit using an asymptotic power law and 
correction term. We find 

In R 2 ~ z sp In t - (t) R r 1/A + C (6) 

with z sp = 0.921(10) and <p R = 1.703. 

It has been argued that 77 = quite generally for extremal dynamics [5,27]. 
Our data for the one-dimensional ABS model support this conclusion: on 
a double-logarithmic plot, n(t) clearly grows more slowly than a power law. 
While 77 = is compatable with n(t) growing without limit as t — > 00 (for 
example, oc (lnt)^, as suggested in Ref. [7]), our results support the conclusion 
that n(t) saturates at a finite value at long times. Specifically, we are unable 
to fit the long-time behavior using an expression of the form n ~ (lnt)^. On 
the other hand, we find d\nn/d\nt oc t~ u , with u ~ 0.149, suggesting that 
n(t) ~ rioo exp[— ct^]. In fact an excellent fit is obtained using c = 1.92 and 
rioo = 14.574, as can be seen in Fig. 6. (Saturation of n(t) does not occur on 
the time scale of the simulation; for the anistropic case, shown in the inset of 
Fig. 6, saturation is in fact evident.) 

In the absorbing phase (q < q c ,E), the survival probability must vanish 
as t — > 00. Our data follow P s ~ t~ 5 e~ t / T , where r ~ \q — q c ,E\~ Un , with 
u\\ = 2.54(2). On the other hand, for q > q cE , the survival probability tends 
to a finite value as t — > 00. We obtain lim^oo P s = P^ ~ (q — q c ,E)^ , with 



10 



P' = 0.20(1). (In DP and allied models (3' = (3 [30], but this need not hold for 
models in other universality classes.) 

In the CP and other nonextremal models, the spread of activity in the 
supercritical regime follows a simple pattern: the size R of the active region 
grows linearly with time, and the number of active sites n grows oc t d . Our 
observation of subdiffusive spreading (and saturation in n) at the critical point 
lead us to investigate supercritical spreading in the ABS model. We find that 
spreading is indeed sublinear. For example, using q = 0.7 in a study extending 
to t ~ 9 x 10 6 to avoid transient effects, we obtain R 2 ~ t x with \ — 1-20(4) 
and n ~ t x with A = 0.61(2). (The exponent governing R 2 should be twice 
that for n, since active regions have a finite activity density in the supercritical 
regime.) Once again, extremal dynamics is seen to limit the growth of activity. 

D. CP3 Model 

We performed spreading simulations of the CP3 model, using the approach 
described in the previous subsection. Each realization is followed up to a 
maximum time of 6 x 10 4 . Using the power-law behavior of P s {t) and n(t) as 
the criterion for criticality, we find q c = 0.63525(3) for the CP3. (Note that 
this is some 5% smaller than the critical value of the corresponding extremal 
model.) Analyzing the local slopes, we obtain 5 = 0.162(2), r\ = 0.312(2) 
and z sp = 1.265(4). These values are fully consistent with those for directed 
percolation (see Table I), confirming that the CP3 model belongs to the same 
universality class as the original contact process. 

A striking difference between extremal and nonextremal models with an 
absorbing state is that the spread of activity in the critical process is much 
slower in the former. This is of course reflected in the value 77 = for extremal 
models, (while for example rj = 0.314 for DP in one dimension), and in the 
subdiffusive growth in R 2 in the ABS model. In Figs. 7 and 8 we compare 
typical evolutions in the ABS model and its nonextremal analog, the CP3, at 
their respective critical points. It is evident that the activity spreads much 
more slowly in the ABS than in the CP3. A further notable difference is that 
in the ABS a site can remain active for a very long time, i.e., while the site 
is not the minimum site or a neighbor of it. Thus the rates of both addition 
and loss of active sites are much smaller in the critical extremal process than 
in the corresponding nonextremal one. 

E. Extremal CP 

In light of the discussion of Sec. II, it is of interest to study the behavior 
of other absorbing-state models under extremal dynamics. As a first step we 
report simulation results for the extremal- absorbing contact process (CPea)- 
We performed spreading simulations to determine q Cj E and the exponents 5, i] 
and z sp , using simulations running to a maximum time of 6 x 10 4 in 5 x 10 5 



11 



independent realizations. We find q c ^E = 0.79415(5) for the extremal CP, 
compared with 0.76733 for the original (non-extremal) process. (Note that, as 
in the comparison between the CP3 and ABS models, g Cj £ > q c . This may 
reflect the slower spread of activity under extremal dynamics.) 

As in the case of the ABS model, the decay of the survival probability at 
the critical point follows an expression of the form of Eq. (5), here with best- 
fit parameters 5 = 0.0855 and <f)p = 0.226. The exponent 5 is essentially the 
same as found for the ABS model, while the correction term is about twice as 
large. At the critical point the derivative d\nn/d\nt ~ t~ 0A , again indicating 
a behavior of the form n(t) ~ exp[— ct'^j, here with uo = 0.1. The mean- 
square distance of active sites from the origin grows in a manner similar to 
that in the ABS model. We are again able to fit the data for R 2 using an 
expression of the form of Eq. (6), with z sp = 0.932 and 0^2 = 2.026. These 
results strongly suggest that the CPea belongs to the same universality class 
as the ABS model. 

We turn now to the rather surprising behavior of the stationary probability 
density p(x) in the extremal CP. Recall that mean-field theory predicts p(x) = 
20 (x — 1/2) for q < q c ^ E = 1/2, while for q > 1/2 there are two steps, one at 
x — k = q 2 / (3q — 1), the other at x — q. In simulations of the CP# on a ring 
we find a single step discontinuity for q < q cE = 0.79415, and, for q > q CtE , 
a pair of steps, one at x — q, the other at x = q s < q. The positions of the 
singularities as obtained in simulation (using data for system sizes L = 100, 
200,..., 1600 to extrapolate the position in the infinite-size limit), are shown 
in Fig. 9. The lower singularity q s is seen to bifurcate from q = q cE just at 
the critical point, in qualitative agreement with MFT. Note however that the 
position of the singularity is not constant for q < q c E , as predicted by MFT. 

The density p(x) is shown for q = 0.794 ~ q c E and q = 0.85 in Fig. 10. In 
the latter case it is evident that the step at x = q is sharp (this is true even 
for small systems) and derives from the singular nature of the updating rule. 
The step at x = q s , by contrast, is subject to finite-size rounding, and becomes 
sharper with increasing system size, as is characteristic of a critical singularity. 
The finite width of the peak at x = q c (in the process with q = q c ,E), appears to 
be a finite size effect as well: it becomes sharper with increasing L, suggesting 
that the singularities merge in the limit L — > oo. 

F. Anisotropic ABS Model 

The scaling behavior of the Bak-Sneppen model changes when the updating 
rule is asymmetric [31]. The same critical exponents are found for a highly 
anisotropic version in which at each step, only the minimal site and its neighbor 
on the right are updated [32], and for weak anisotropy [16,33], so that one may 
identify an anisotropic BS universality class. In this section we report results of 
spreading simulations of the anisotropic absorbing BS model. To obtain these 
results, we simulated the anisotropic ABS (in the highly anisotropic version) 
in studies extending to a maximum time of 1.6xl0 5 , using 3xl0 5 realizations. 
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Analysing the local slope S(t) = d In P/d In t, we determined the threshold 
of this model to be q Ct E = 0.72370(2). (The best previous estimate is 0.7240(1) 
[34].) A typical evolution of the critical spreading process is shown in Fig. 
11. The local slope 5(t) yields the estimate 5 = 0.234(5). For anisotropic 
models we define R 2 {t) as the mean-square radius of gyration, i.e., the distance 
is measured relative to the current center of mass of the set of active sites, 
rather than to a fixed origin. This is done to eliminate a spurious contribution 
due simply to the overall drift in the active region. For the anisotropic ABS 
model R 2 {t) may again be fit with an expression of the form of Eq. (6), with 
z sp = 1.425(10) and 4> R = 2.3(2). The exponents 5 and z sp are quite different 
from those of the isotropic model. Despite these differences, we again find 77 = 
for the anisotropic model. As before, the mean number of active sites n(t) 
saturates at long times, more rapidly in fact than in the isotropic ABS model 
(see Fig. 6, inset). We are able to fit the data well using n(t) = n^l — e _rfl/4 ] 
with parameters = 5.206(3) and c = 0.348. 

The nonextremal model corresponding to the anisotropic ABS model is 
a two-site contact process, CP2, which is simply the CP3 with updating re- 
stricted to the central site and its neighbor on the right. We have verified that 
the spreading exponents of the CP2 model are those of directed percolation. 
(Here again, we define R 2 as the mean-square radius of gyration.) This leads to 
the interesting conclusion that a perturbation (asymmetric updating) that is 
irrelevant for a nonextremal model is relevant for the corresponding extremal 
system. (We note that, because the two sites in the CP2 are updated in the 
same manner, the model does not fall in the so-called anisotropic-DP class, for 
which bonds along different axes are present with different probabilities [35].) 

IV. CONCLUSIONS 

We investigate the relation between extremal dynamics, exemplified by the 
Bak-Sneppen model, and nonextremal models exhibiting a phase transition to 
an absorbing state, using general arguments, mean-field theory and simula- 
tion. The relation between the BS model and directed percolation was already 
suggested some time ago [4,5]. Here we clarify this connection by showing 
how a generic absorbing-state model can be transformed to an extremal one 
via the associated extremal-absorbing model. The nonextremal 'precursor' of 
the BS model is a three-site contact process [24], CP3, which, like the original 
CP, belongs to the directed percolation universality class. The BS model and 
the extremal version of the CP belong to a common universality class that 
may be dubbed 'extremal DP' (EDP). A number of extremal dynamics classes 
distinct from EDP are discussed in [5]; another example is the anisotropic BS 
model. We expect that further extremal dynamics universality classes exist, 
for example an extremal parity-conserving class [36], although examples of the 
latter have yet to be studied. 

Our results for the critical exponents of the EDP class, which includes the 
BS model and the extremal CP, are compared against those of ordinary DP (in 
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one spatial dimension) in Table I. (Here we have taken the values 77 = and 
/3 — 1 to be exact for EDP.) The differences between the two sets of exponent 
values are evident. Our results r = 1 + 5 = 1.084(1), and z = v\\/v±_ = 
2.12(1) are in agreement with the earlier estimates [5] of 1.07(1) and 2.10(5), 
respectively. Our result is however somewhat higher than Grassberger's result 
r = 1.073(3) [7]. 

Certain scaling relations are expected to hold among the critical exponents 
[21,23,30]. In spreading processes one expects z sp = 2v>±/i>\\; our data are 
nearly consistent with this, yielding 2v±/i>\\ — z sp = 0.022(14). The relation 
P' = 5v\\ is also satisfied: our data yield f3' — 5v\\ = —0.013(14). Finally, we 
consider the generalized hyperscaling relation [37] 

2 { 1 + J>) 6 + 2r] = dZsp ' (7) 

in d dimensions. Using our data, we find the difference between the two sides 
of this relation to be 0.09(6). Our results are marred by another inconsis- 
tency that may reflect corrections to scaling or finite size effects: the product 
(P/ l/ ±)~ 1 ( l/ \\/ l/ ±) l '\\ 1 : with the first two factors determined from finite-size scal- 
ing at the critical point, and the final factor obtained from the decay of the 
survival probability in the subcritical regime, should equal (5 — 1; our data 
yield 1.08(3). These minor inconsistencies suggest that one or more of the 
exponents may be in error by 5% or so. Refining their values will require accu- 
mulating larger data sets in simulations of larger systems, a task we leave for 
future work. (The studies reported here were quite demanding computation- 
ally, representing approximately 6 months' cpu time on an alpha workstation.) 

In the course of our study we revisit a three-site contact process (CP3) 
that is the nonextremal analog of the BS model [24]. We verify that the 
CP3 belongs to the universality class of directed percolation, as expected [7]. 
Similarly, we define extremal and extremal-absorbing versions of the original 
contact process (CP^ and CP ea, respectively) and verify that their scaling 
properties are the same as those of the BS model. The stationary probability 
density for the CP^ follows, in general terms, the predictions of mean- field 
theory, but certain interesting differences exist, as detailed in Sec. III.E. 

It is clear that when an absorbing-state model is modified to follow extremal 
dynamics, its critical exponents are altered. In general extremal dynamics tens 
to slow the spread of activity in the critical and supercritical regimes. One may 
nevertheless inquire whether any more general features of the original model 
are preserved under 'extremalization'. A candidate for such a conserved prop- 
erty is the critical dimension d c . In critical phenomena, various universality 
classes (differing in the symmetry group of the order parameter, or the pres- 
ence of conserved quantities) may share the same d c if the algebraic structure 
of their continuum description (in particular, the power of the lowest-order 
nonlinear term in the order parameter, in a Landau- Ginzburg- Wilson effective 
Hamiltonian) is the same. Thus d c — 4 for all models in the n-vector family. 
Extending this idea to extremal models is questionable, since there is no con- 
tinuum description at hand. (At first glance, the notion of extremal dynamics 
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in a description using a continuous activity density seems problematic, since 
there is always one and only one extremal site.) Be that as it may, it seems 
plausible that if the field theory for DP [18,19,38] were somehow modified to re- 
flect extremal dynamics, the dominant nonlinearity would not change, so that 
d c would retain its value of four, as in DP. The upper critical dimension d c = 4 
for the BS model was established some time ago by Boettcher and Paczuski 
[39]. Our argument suggests that extremal versions of other absorbing-state 
models have the same upper critical dimension as the corresponding nonex- 
tremal model. We hope to test this prediction in future work. 

Studying the anisotropic ABS model and its nonextremal counterpart, the 
CP2 model, we find that anisotropy is a relevant perturbation for extremal 
DP, while it is irrelevant for the corresponding nonextremal class. We suspect 
that other perturbations, such as diffusion, may exhibit a similar pattern of 
relevance. 
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APPENDIX: MEAN-FIELD THEORY 



A. Extremal dynamics as a zero-temperature limit 

There are several ways of formulating a mean-field theory (MFT) for ex- 
tremal models. First we consider an approach [11,33] in which the probability 
of a site i being chosen as the central site is proportional to e~^ Xl \ extremal 
dynamics is recovered in the limit (5 — > oo. (In the present discussion the 
parameter (3 bears no relation to the critical exponent denoted by the same 
symbol in the main text.) Applied to the BS model, this approach yields the 
distribution p(x) = (3/2)Q(x - 1/3)0(1 - x) when (3 -> oo [11,33]. 

In the ABS model the distribution p(x) evolves via 



= - e -^p{x)Q{q -x)+3j q e-f»p(y, t)dy - 2p(x, t) j" e^Mv, t)dy 

(8) 

The first term represents a site with value x being selected as the central site, 
which is only possible if x < q. The second term reflects updating three sites 
with new variables uniform on [0, 1], with the integral representing the overall 
rate of events. The final term represents updating of the two neighboring sites, 
and is obtained using the mean-field factorization of the nearest-neighbor joint 
probability density: p(x,y,t) ~ p(x,t)p(y,t). Note that in writing Eq. (8) we 
have associated a time increment dt = 1/N, with N the number of sites, with 
each event. 

Eq. (8) admits an infinite set of stationary solutions for which p(x) = on 
the interval < x < q. These represent absorbing states. To seek an active 
stationary solution we let / = e~^ y p(y,t)dy, and equate the time derivative 
to zero, yielding 

= 2Z + e(g-*)e-g' (9) 

To determine I we multiply Eq. (9) by e~^ x and integrate from x = to x = q, 
leading to / = (e^/ 3 - e~^)/[2(l - e"^ 3 )], so that 

3/>-/3/3 _ p -/3q\ 

v( x ) = V? - I MO) 

P{) (e-^-e-^) + e(q-x)e-^(l-e-^) 1 } 

In the limit (5 — > oo, we find, for q > 1/3, the singular density p(x) = 
(3/2)0 (a; — 1/3). This is precisely the MF result for the original BS model. 
[When we take f3 — > oo, the above expression reduces to p(x) = (3/2)0(x — q) 
for q < 1/3. But this density is not normalized on [0, 1] and so must be re- 
jected. We are left with only absorbing stationary solutions for q < 1/3.] Thus 
Qc,e = 1/3 in the MFT of the absorbing Bak-Sneppen model. Note that the 
parameter q is irrelevant for q > q c E = 1/3, as was argued in Sec. II. 
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A moment's reflection shows that the evolution of p(x) in the (nonextremal) 
CP3 model is given by Eq. (8) with /3 — 0, since all active sites are then equally 
likely to be chosen as the central site. Taking the limit j3 — > of the stationary 
solution, Eq. (10), one finds, for q > 1/3, the stationary density 



p(x) 



|(3 - q x ), x < q 



3 
2> 



q < x < 1 



Eq. (11) confirms that the stationary density of a (nonextremal) model ex- 
hibiting an absorbing state phase transition is characterized by a steplike sin- 
gularity, as asserted in Sec. II. For q < 1/3, Eq. (11) yields an unphysical, 
negative density, showing that q c = 1/3 for the CP3, in the mean-field approx- 
imation. 

The foregoing analysis is readily extended to the extremal-absorbing con- 
tact process (CP ea) defined in Sec. II. The rate of events is again given by 
/ = Jq e~ l3x p(x)dx. At each event, there is a probability of 1/2 that the 
central site (which must have x < q) is replaced, while with probability 
1/2 a neighbor is updated. Thus the loss terms in the equation for p(x,t) 
are: — (l/2)[e~ l3x Q(q — x) + I]p(x). The gain term corresponding to up- 
dating of the central site is simply 1/2, but for updating a neighbor it is 
(J/2)[l — P(q) + (l/q)Q(q — x)P(q)], where P(x) = J* p(y)dy is the probabil- 
ity that a given site i has Xi < x. (The reason is that when updating an active 
neighbor, the new variable is chosen from the distribution uniform on [0, q].) 
Thus the MF equation of motion is 



dp 
di 



1 



l3x p(x)Q(q - x ) + - 



2 - P(q) + 



Q(q — x) 



P(q) -p(x) 



12) 



To find the stationary solution p(x) we first note that for < x < q, setting 
dp/dt to zero yields 



p(x)=AI/(I + e- fix ) 



(13) 



where A = 2 + (q 1 — l)P(q). Integrating Eq. (13) from x = to x = q, we 
find 



„ / s . A f 1 du 
P(q) = Aq-- - 



+ u 



(14) 



If we now multiply Eq. (13) by e @ x and integrate over the same interval, we 
obtain 



A f 1 du 



(3 J e -f3q I + U 



(15) 



leading to P(q) = 2 — q 1 in the stationary state. We see that g Cj £ = 1/2 as in 
the MFT of the original contact process. Evaluating the integral in Eq. (15) 
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one finds I = (e f3K - e /3g )/(l - e /3k ), where re = q 2 /(3q - 1). The stationary 
density is 



p{x) 



1 (e -^-e^)[l + (2-g' 1 )9(g-x)] 
q (e~P K - e~fo) + e"^(l - e-P K )Q(q - x) 

For q > 1/2, we have re < q, and in the limit (3 — > oo, 

'0, X < K 

p{x) = 



(16) 



i 



, k < x < q 
x > q 



(17) 



which is normalized and exhibits step-function singularities at x = re and 
x = q. For g < 1/2 on the other hand, re > q and Eq.(16) does not yield an 
acceptable probability density, and we conclude that the only stationary state 
is the absorbing one. The critical point of the CPea thus falls at q Cj E = 1/2 
in MF approximation. 

Taking (3 — > in Eq. (16), we obtain the probability density for the original 
CP: 



p(x) 



2q-l 



, x < q 



; , x > q 

Finally, for the extremal CP, the equation of motion is 



(18) 



dt 2 P[ ] + 2 



2-P(q) + ^-^-P(q)-p(x) 



(19) 



with / = e /3x p(x)dx. To find the stationary solution we write 

_ 2 + [q^ejq - x) - l]P(q) 
P[X) ~ 1 + J-i e -/fc 

Integrating from to q and solving for P(q) we find 

2(<jr- 7 ) 



9 + 7(9 1 - 1) 



where 



1, /+ 1 
7 = — In 



/3 / + e-Pi 

Now multiply Eq. (20) by e~@ x and integrate from to 1 to obtain 



(20) 



(21) 



(22) 



1 = A 



e @ x dx 
I + e-P* 



+ A' 



, f 1 e P x dx 
I + e-P* 



(23) 
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where A = 2/[q + 7(g~ 1 — 1)] and A' = 'jA/q. If P(q) > 0, the first term on 
the r.h.s. of Eq. (23) is nonzero and dominates as (3 — > oo. Equating the first 
term to unity then leads to 7 = q 2 / (3q — 1) = k, and thence to P(q) = 2 — q^ 1 
which is positive for q > 1/2. A simple calculation then yields the distribution 
of Eq. (17) in the limit f3 -> 00. 

If q < 1/2 the above solution is not valid since it implies P(q) < 0. We 
therefore take P(q) = 0, implying 7 = q, and so A = A' = 2. Eq. (23) now 
reads 

9 74- e ~ l3q 

Solving for / and inserting the result in Eq. (20), we find in this case 
lim^oo p(x) = 2Q(x — 1/2). These results have been verified via numerical 
integration. 



B. Extremal dynamics on a complete graph 

Another approach to formulating MFT for the BS model considers extremal 
dynamics on an iV-site complete graph or random-neighbor model (two neigh- 
bors are selected at random each time a site is updated); the stationary density 
p(x) becomes a step function in the infinite-size limit [12,16,40,41]. We now 
extend this approach to the ABS model. Let P(x) = Probfxj < x] = J* p(y)dy 
be the distribution function and let Q(x) = 1 — P(x). By definition Q is a 
non-increasing function with Q(0) = 1 and Q(l) = 0, since p(x) = outside 
the interval [0,1]. 

Activity in the ABS model is predicated on the minimal site x m i n being 
smaller than q; the probability of this event, under the MF factorization, is 
1 — [Q(q)] N . Given x min < q, updating the extremal site and two neighbors 
results, on the average, in the increment: dP(x) = (1/N){— [1 — Q(^<) JV ] — 
2P(x) + 3x}, where a;< = min{x, q}, so that the first term represents loss of 
the minimal site, the second removal of two neighbors, and the third random 
replacement of the three site variables with numbers uniform on [0,1]. If we 
adopt a time increment dt = 1/N for each such event, the equation of motion 
for P is 



dt 



= _[i _ Q( x< , t ) N ] + [1 - Q(q, ty][3x - 2P(x, t)] (25) 



Note that the evolution ceases if Q(q,t) = 1, i.e., if there are no active sites. 
(Since Q is nonincreasing Q(q) = 1 =>- Q(a;<) = 1.) 

For q > 1/3, the stationary solution to Eq. (25) corresponds to a density 
p(x) that approaches a step function, (3/2)Q(x — 1/3), as iV — > 00. A simple 
calculation yields the dominant contribution for large iV: 

Q~{ {1 - 3X)1/N > t X<1 * (26) 

V ~ \ |(1 - x) + 0(e- const - N ) x > § 1 ' 
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(One should note however that the convergence is nonuniform in x, being 
slower the closer x is to the critical value of 1/3.) For q < 1/3 we are unable 
to find an accepable stationary solution with Q < 1 (i.e., p > 0), for x < q, 
and conclude that only absorbing solutions exist. 

The analysis of the ABS model on a complete graph confirms that in the 
infinite-size limit, the model enjoys the usual properties of the BS model for 
q > Qc,e = 1/3, and falls into the absorbing state for q < 1/3. 

The evolution of P(x,t) in the extremal CP follows, in MF approximation, 
the equation 

= - Q(x,t) N ] + |- Ip(x,t) + \[x*P(q,t)+xQ{p,t)] (27) 

where x* = mm{x/q, 1}. Numerical integration shows that the solution con- 
verges, for large N, to a stationary distribution consistent with the singular 
density found above in the limit f3 — > oo. 
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TABLES 

TABLE I. Critical exponents for the one-dimensional absorbing Bak-Sneppen 
model (ABS) and contact process (CP). CP exponents from Refs. [23] and [41]. 



Exponent 


ABS 


CP 


P 


1 


0.27649(4) 




0.20(1) 


(=P) ' 


»\\ 


2.54(2) 


1.73383(3) 


P/"± 


0.77(1) 


0.25208(5) 


v \\/ v ± 


2.12(1) 


1.58071(11) 


5 


0.084(1) 


0.15947(3) 


7] 





0.31368(4) 




0.921(10) 


1.26523(3) 



TABLE II. Spreading exponents for the CP3 and CP ea models and the 
anisotropic absorbing Bak-Sneppen (AABS) model in one dimension. 



Exponent 


CP3 


CPea 


AABS 


5 


0.162(2) 


0.0855(20) 


0.234(5) 


7] 


0.312(2) 








Z.sp 


1.265(4) 


0.932(20) 


1.425(10) 
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FIGURE CAPTIONS 



FIG. 1. Stationary activity density (filled symbols) p versus system size L in 
the one-dimensional ABS model at the critical point. Open symbols: difference 
between p and the fitting function, Eq. (3), shifted vertically for visibility. 

FIG. 2. Mean lifetime r (filled symbols) versus system size L in the one- 
dimensional ABS model at the critical point. Open symbols: t/L u ^^ u± (shifted 
vertically for visibility). 

FIG. 3. Moment ratio m for the ABS model versus system size L~ a25 . Points: 
simulation data; line: best linear fit, m = 1.0295 + 0.268L~ 25 . 

FIG. 4. Local slope S(t) versus \jt in the ABS model, q values (bottom 
to top) 0.66700, 0.66701, 0.66702 and 0.66703. Inset: data for q = 0.66701 
plotted versus 1/t - 25 . 

FIG. 5. Survival probability P s (t) in the ABS model at the critical point, 
q = 0.66701. The nearly constant function represents the ratio of P s to the 
fitting function, Eq. (5). 

FIG. 6. Mean number of active sites n(t) in the ABS model at the critical 
point, q = 0.66701. The solid curve represents the fitting function described 
in the text. Inset: a similar plot, for the critical anisotropic ABS model. 

FIG. 7. Spread of activity in a typical realization of the critical CP3 model 
(q = 0.63525). 

FIG. 8. Spread of activity in a typical realization of the critical ABS model 
(q = 0.66701). 

FIG. 9. Position q s of the singularity in the stationary probability density of 
the extremal CP. The upper line of singularities, x — q, bifurcates from q s at 
the critical value q Ci E- 

FIG. 10. Stationary probability density p(x) in the CP# for q = 0.794 (left) 
and q = 0.85 (right); system size L = 1600. 

FIG. 11. Spread of activity in a typical realization of the critical anisotropic 
ABS model (q = 0.72370). 



24 




500 



X 



0.084 ■ r 

-0.085 



0.085 



0.086 



0.087 



0.088 




0.0000 0.0001 0.0002 0.0003 0.0004 

t" 1 



0.0 



-0.2 



-0.4 



-0.6 



-0.8 



■1.0 



i 1 i 1 r 



~i 1 r 



i 1 r 




-1.2 



J i I i L 



J I I I I I L_ 



J i L 



2 4 6 8 10 12 14 



8 



Q. 4 



t i i r 




_1 I I L_ 



111 

0.70 0.75 0.80 0.85 0.90 




0.6 0.7 0.8 0.9 1.0 




4 5 6 7 8 9 10 11 

In L 



20 
18 
16 
14 

12 - 

10 - 

8 



